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ABSTRACT 


The optimal control of 
maneuvers and vibrations of 


three-dimensional large-angle rapid 
a Shuttle-mast-reflector system is 


considered. The nonlinear equaions of motion are rormuiated by 
using Lagrange's formula, with the mast modeled as a continuous 
beam subject to three dimensional deformations. The nonlm-ai 
terms in the equations come from the coupling between the angular 
velocities, rhe modal coordinates, and the modal rates. 

Ponnry&gin’ s Maximum Principle is applied to the slewing problem, 
to derive the necessary conditions for the optimal controls, 
which are bounded by given saturation levels. The resulting two- 
point boundary-value problem is then solved by using the 
quasi linearization algorithm and the method of paiticular 
solutions. The numerical results for both the linearized system 
and the nonlinear system are presented to compare the differences 
in their time responses. 

The study of the large-angle maneuvering of the Shuttle- 
beam-reflector spacecraft in the plane of a circular earth oibit 
is extended to consider the effects of the structural offset 
connection, the axial shortening, and the gravitational torque on 
the slewing motion. The offset effect is analyzed by changing 
the attachment point of the reflector to the beam. As the 
attachment point is moved away from the mass center of the 
reflector, the responses of the nonlinear system deviate from 
those of the linearized system. The axial geometric shortening 
effect induced by the deformation of the beam contributes to the 
system equations through second order terms in the modal 
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amplitudes and rates. The gravitational torque effect is 
relatively small. 

Finally the effect of additional design parameters {such as 
related to addit . :>nai payload requirement) on the LQR based 
Qf orbiting concrc i/stmctural system is exa*>ii.ied. 
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optimal multicriteria 


r a : n t s not only include the limited mast 
“he system but also in.lude the control 
is transiert response time of the system) . 
are derived for minimizing the cost 


function and setting the variation of the cost function witn 
respect to the design variables to zero. The simple models ■ - 
uniform solid and tubular beams are demonstrated here with two 
typical additional payload masses: (1) symmetrically distrirut- 

with respect to the center of the beam; ; 2) asymmetrically 
distributed with respect to the center of the beam. For the 
s 1 id and tubular beams, the length and material properties are 
assumed equal. By considering the transient response of pitch 
angle and free-free beam deformations in the orbital plane, the 
optimal outer diameter of the beam and all feedback control can 
be determined by numerical analysis with this multicriteria 
approach . 
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I. INTRODUCTION 


The present grant, NSG-1414, Suppl . 12 extends the research 

initiated in May 1977 and reported in Refs. 1-16' for the grant 
years May 19 7 7 - May 1999. This effort has focused on the 
problem of shape and orientation control of large, inherently 
flexible proposed space systems. Possible applications proposer 
f c r these large space systems ( LSS ) include: Earth, observation 

and resource sensing systems; orbitaily based electronic mai^. 
tra:ism_ssion; large scale multi-beam antenna communication 
systems; as platforms for orbital-based telescope systems; and as 
m-orbit test models designed to compare the performance of LSS 
systems with that predicted based cn scale model Earth-based 
1 aoor story experiments and/or computer simulations. In the last 
several years the grant research has focused on the orbital model 
of the Spacecraft Control Laboratory Experiment (SCOLE).' 

The present report is divided into five chapters. Chapter 
II is based on a paper presented at the 1990 International 
Conference on the Dynamics of Flexible Structures in Space and 
describes rapid three dimensional maneuvers and vibration 
suppression of the asymmetrical flexible SCOLE configuration. 
Pontryagin's maximum principle is applied to both the linearized 
and nonlinear system equations to develop the necessary 
conditions for the optimal multi-control problem. The resulting 
two point boundary value problem is then solved based on the 

’ References cited in this report are listed separately at the end 
of each chapter. 
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quasi linear! z at ion algorithm and the method of particular 
solutions. The numerical results for both the linear and 
nonlinear systems are presented to compare the differences m 
their tine responses. 

Chapter III is motivated by a paper presented at the 17th 
International Symposium on Space Technology and Science m which 
the effects of the structural offset for asymmetrical 
configurations (sues, as SCCLE ; , axial shortening, and 
gravitational torque during a rapid slew are evaluated. For this 
study a two dimensional model of the SCOLE Shuttle - (flexible; 
be am - r e f 1 e c tor system is considered. Both lineai and nonlinear 
system models are treated. 

In Chapter IV the effect of additional system design 
parameters (such as those related to the placement of additional 
payloads) on the LQR based design of an orbiting control/ 
structural system is analyzed. This multicriteria numerical 
optimization approach is considered for minimizing an LQR type 
cost function where the system design parameter is the outside 
diameter of a solid and/or tubular beam, subject to constraints 
on the total system mass, control saturation levels, and 
transient settling time. Different combinations of additional 
payload masses are considered. 

Finally, Chapter V describes the main general conclusions 
together with general recommendations. The thrust of this effort 
has been redirected to provide more direct support to the new 
NASA Controis/Structures Interaction Program (CSI), particularly 
as evidenced by Chapter IV, and our follow-on proposal, Ref. 18. 
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II. 


OPTIMAL LARGE ANGLE MANEUVERS OF A FLEXIBLE SPACECRAFT 


1. INTRODUCTION 

Many authors have considered the problem of large-angle rapid maneuvers 
of flexible spacecraft [1-5]. The direct solution of the open-loop two-point 
boundary-value problem (TPBVP) for three-dimensional (3-D) slews of flexible 
spacecraft resulted in numerical problems with rank.— deficient matrices as 
stated by Chun [3]. However, a different numerical method may be used to 
overcome this difficulty. In this paper, the problem has been solved 
successfully by using the quasilinearization algorithm and the method of 
particular solutions for 3-D slews of an asymmetrical flexible spacecraft, 
namely, the Spacecraft Control Labratory Experiment (SCOLE) configuration. 

The open-loop slewing approach has several obvious distinct properties. 
First, the control law is easy to implement in practice for both ground tests 
and space flight tests. Second, the open-loop solution may serve as a good 
reference for the feedback control law design, as proposed by Chun [3], and 
Meirovitch [4], in which the open-loop solution for a rigid (instead of a 
flexible) spacecraft is used as the nominal reference trajectory. As an 
extension to Refs. [3] and [4], it may be helpful if the open-loop solution 
for the 3-D slew of a flexible spacecraft system could also be used as a 
nominal reference solution. In addition, through the present study, we can 
also see how different are the responses of the nonlinear system from those of 
the linearized system, and the differences between the flexible and rigidized 
systems . 


2. FORMULATION OF THE STATE EQUATIONS 
2.1 Sy stem Conf iguration 

As shown in Fig. 1, the orbiting SCOLE configuration [6] is composed of a 
Shuttle, a flexible mast, and a reflector antenna. Both the Shuttle and the 
reflector are considered to be rigid bodies. One end of the mast is fixed to 
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the Shuttle at its mass center, o , while the other end is firmly connected to 
the reflector at an offset point, a^ . 

Three coordinate systems, (i j k ), (i j k ) , and (i j k ) , 

OOO 8 8 8 ii» 

representing the orbit's local vertical/horizontal reference system, the 
Shuttle body axis coordinates, and the reflector axes, respectively, are 
adopted in Figure 1. The mass center of the reflector, O f , is located at (x^ 
y ) in the reflector axis system. Three Euler angles (0^ , 0 2 , 0 3 ) or four 

r 

quaternions (q Q q^ q z q 3 ) are used to describe the attitude of the Shuttle 

with respect to the orbiting reference system. 

The undeformed mast is assumed to be oriented along the z axis of the 
Shuttle coordinate system. The 3 _ D deformation of the mast consists of two 
bending deflections U(z,t) and V(z,t) in the x-z and y-z planes, respectively, 
and torsion <f>(z,t) about the z axis. It is assumed that these deformations are 
small as compared with the length of the mast and can be expressed by the 
following modal superposition formula [7]: 

U=E£. (z)a i (t), V=p 7 i (z)a. (t), <fr=EC. (z)^ (t) , (1) 

where £ n . and C are modal shape function vector components normalized by 
a common factor, and a is a scaled modal amplitude associated with the ith 
mode. The free vibration of this structure can be considered as a space 

free-free beam vibration problem with boundary conditions including the masses 
and moments of inertia of the Shuttle and the reflector. The partial 
differential equation formulation for this problem [6,7] can be solved by 
using the separation of variables method. The first five natural frequencies 
and mode shapes have been obtained by Robertson [7], and will be used in this 
paper . 

2.2 Kinetic Energy 

The kinetic energy of the system about the mass center of the system, c, 
can be expressed as 
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( 2 ) 


T = U i;i 2 d» ♦ |[|j;i 2 dm - 4r(lj ; l do ) Z ] 

* i[J r l'l 2<ta - -s;(J r I^H 2 ] ' -sr(J b H-(I r H 

_ T + x + T - T 

s b r c 

where r is the position vector from o^ to an arbitrary mass element in the 
system and m is the total mass of the system. The integration subscripts, 
"s", "b” , and "r", mean that the corresponding integration is throughout the 

Shuttle, the beam, and the reflector, respectively. 

Kinetic Energy of the Shuttle The first term in Equation (2) is the 
rotational kinetic energy of the Shuttle about o^ , 

T = -^i) T J <i) (3) 

9 2 s 

where <d is the matrix describing the angular velocity vector of the Shuttle, 

<0, and J is the inertia matrix of the Shuttle. 

8 

Kinetic Energy of the Mast As Shown in Figure 2, the position vector of an 
element dm and its velocity are, respectively, 

r=b+p , b=Ui +Vj +zk (4) 

8 8 8 

r=v + G) (z)x p (5) 

b b 

where p is a vector within the cross section of the beam and 

v =Ux +Vj + G) x b (6) 

b 8 8 

G) (z) is the angular velocity of the element , 
b 

G> (z) = (3+<|>i+<i>j+^k (7) 

b x 8 y s z 8 

where =-(dV/dz), <p =6 U/dz , and 

x y z 
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During the integration process of equation (2), one needs to do the 
following calculations, 


J |u b xp| 2 dm=j |p« b [j (p T pE-pp T )dAjo> b |dz 


( 8 ) 


where <0 is the matrix representation of o> and E is an identity matrix. For 
b b 

the circular cross section of the beam assumed here, the inner integration in 
Equation (8) turns out, in the local beam coordinates , to be. 


J (p T pE - pp T )dA.j - 1^ 


-00 

2 

0-0 

2 

0 0 1 


= J (r) 


(9) 


where I is the polar moment of inertia of the beam. Due to the small 
P 

deformation of the beam, the local beam axes is assumed to be related to the 
Shuttle axes by a transformation matrix, 


1 <t> 

z y 

<f> 1 ^0 

Z X 

$ 1 
y * 


( 10 ) 


Therefore, J (r) in Equation (9) can be transformed to the Shuttle axes by the 
following similarity transformation: 


.<•> 


. R" r J (r> (R* r ) T 


(ID 


After substituting all related terms into the second term of Equation (2), and 
neglecting all the third and higher order terms in the modal amplitude vector, 
a, the modal rate vector, a, and their coupling, one can arrive at 


T =\) T J u + -ia T I a +u T h 
b 2 b 2 b b 


( 12 ) 


where 1 is a constant matrix. The elements of matrix J and the vector h 
b b b 

have the following forms. 


(J ) =c *ct T m +b T M a, (h ) =ct T g +o T G a 
b ij ij ij ij b'i i 


(13) 
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where , m^, * ^ij' an< ^ are constants * 

matrices , respectively * 


constant vectors , and constant 


Kinetic Energy of the Reflector After using a development process similar to 
that for the mast, one can obtain, 

T =\} T J 0 ) + -!(x T I a *ti) T h (14) 

r 2 r 2 r r 

Here the inertia matrix of the reflector needs to be considered in the 
development process . 

The Coupling Term in Equation (2) can also be written in the form of Equation 
(12) for consistence, 


T =\) T J <i ) + — ct T I a +(*) T h 

c 2 c 2 c c 


(15) 


After substituting Equations (12,14,15) into Equation (2), one obtains. 


T* { J +J +J } it) + “« T I oc */{ h +h } 

2 0 1 2 2 1 Z 


(16) 


where J , J , J are 3x3 matrices, h , h are 3x1 vectors, I is an nxn 
0 12 12 

constant matrix, and 


( J o ) inconstant, ( Vi j® 


(h ) = d T g , (h ) = d T G a 

v i ' i *i ' 2 i i 


(17) 


where ■ , , M^, and G t are constant vectors and matrices. 


Potential Energy 

The elastic potential energy of the flexible part is 


8 *} 


(18) 


where El is the bending stiffness and G is the modulus of rigidity of the 
beam. 
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Generalized Forces 

The virtual work down by the controls is 

4 

5 W=N .3e+ E 5 • & r t (19) 

l j-2 j J 

where 30 is a variation vector which has the direction of an axis of rotation 

of the Shuttle and the magnitude of the rotation angle about this axis, Sr^ 

are the virtual displacements of the position vector r at the location of the 
controls, N is the control torque vector on the Shuttle, and are the 
control force vectors on the beam and the reflector. The shift of the center 
of mass of the system is also considered in this development. After expanding 
these terms by using the associated relations and the transformation matrix 
(10), dropping all third and higher order terms in a, one can get, 

S IMS’lVjkVjl* <20) 

where f =[f f t f t ] T , f =[ f ^ f^ ] T , are control variables; the elements of 
matrices T and T have the similar form as those in Equation (17), and up to 

the first order terms in a have been retained for later use. Then, the 

generalized forces are, 


Q =f + E T f 
0 l j-2 J J 


Q = E ¥ f 
a J-2 j J 


Dynamical Equations 

After constructing the Lagrangian, L=T-V, and using 


d 

It 




( 21 ) 


( 22 ) 


one can get the following equations by discarding second order terms in a, a, 
and their coupling, and retaining a constant mass matrix represented by the 
coefficients for <5 and a, 


J « +[g g- g ] T a =[m*M(a)]S ♦H(a)« +<k (23) 

0 12 3 ” 

[gj g 2 g 3 l w ♦ la =[n+N(o)]S +G(a)« + Q a “*o (24) 
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where 


"=[< w z «3 "z " 3 “3", ".“3] 


(25) 


After finding the inverse of the mass matrix, a state form of Equations 
(23) and (24) can be obtained as. 


u 


0 


= (A +B )5 +(C 0 )u +Doc +(E +F )u 

a p ct 


F g a] , with A, B j , C ^ 


D, E, F^ being constant 


(26) 


where 0=a, B -[B a | B a }•••{ B a], C_-[C 8 ] C B j C 0], and F -[F o i F a 

0C 1 2 6 p 1 Z J w 1 * 


trices; and u-[f f 

lx ly 

f j£ f If f 1 T For the purpose of comparison, the dynamical 
2y* 3x 3y ' 4x 4y 

equations for the rigidized (rigid) spacecraft can be obtained by deleting all 
terms related to <x and 0, this is, 


f If 
1 z * 2x 


<*> 


A o> + E u 


(27) 


where A and E are 3x6 and 3x9 constant matrices, respectively. A linearized 
form of Equation (27) can also be obtained by deleting all nonlinear terms, 


y = Dot + Eu 

The kinematic equations for the quaternions are 


(28) 


1 -V 

q -_«q. 


where 


u = 


0 


^2 

^3 

<*) 

l 

0 

W 3 

^2 

W 2 

^3 

0 

“l 

<*) 

Cl) 

-V). 

0 


(29) 


where q is the 4x1 quaternion vector. 
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3. DERIVATION OP THE OPTIMAL CONTROL PROBLEM 


Necessary Conditions 

A quadratic cost function is used, 

1 r 1 

J=i[ f (a T Q a -*t) T Q w +0 T Q 0 +u T Ru)dt (30) 

4J Q 1 2 3 

where # and R are weighting matrices, t f is the given slewing time. The 
magnitudes of the controls, u, are bounded, 

l u i I* u ib' i=1 »* * • - 9 - od 

The following technique is used to solve this problem [ 8 ]. Pirst, the 
necessary conditions based on Equations (26, 29, 30) are derived. Then, the 
constraints (31) are imposed on these necessary conditions to modify the 
controls • 

The Hamiltonian of the system is, 

H=^(a T Q a +a> T Q <«> +0 T Q B +u T Ru +p T 5q) +y J P 

♦X T [(A +1 )S +(C_)w +Dot +(E +F )u] (32) 

ot p oc 

where p, 7 , and X^ ] T are the costate vectors associated with q, a, o>, 

and p, respectively. By using the Maximum Principle, the necessary conditions 
for the unrestricted optimal control problem are the dynamical equations (26, 
29) plus the following differential equations for the costates. 


• _ dH _ 1 ~ 

P ‘"aq ~ 2 - P 

(33) 

* =-§a = *V -* T * 

(34) 

/"IS = “ Q 2 w "i [qlp -[* T <A+i a )]« -(Cp) T x 

(35) 

V"|? = "V 

(36) 
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and the optimal control. 


I? =0, ==> u=-R‘ 1 (E+F ) T X (37) 

ou a 

The control resulting from Equation (37) is then modified by the 
following saturation considerations. 


u 

i 


-u 

ib 


u 

lb 


if u £ 
i c 



if u i 

ic 



otherwise, 

U =u =-[R _1 (E+F ) T Xl 
i ic Of i 

i-1, 9. 


(38) 


By substituting the control expressions into Equations (26) and (34), one can 
obtain a set of 4(n+3)+2 differential equations for the states and the 
costates. To obtain the control, u, one needs to solve these equations with 
the given conditions: q(0), a(0), w(0), a(0), and q(t f ), a(t f ), <*)(t f ), a(t f ). 


4. Two-Point Boundary-Value Problem (TPBVP) 


One way of obtaining the optimal control law is to transform the above 
necessary conditions into the following TPBVP. Let x represent the state 
vector, and X represent the costate vector. After substituting the control 
expressions (38) into equations (26) and (34), one can obtain two sets of 
ordinary differential equations for the states and the costates, 

(39a) 
(39b) 

with the following boundary conditions, 


''V*’ X) <7.2,,>„1 


i * 2 (*. x ) (7 , 2n) „, 


x(0) and x(t f ) prescribed, X(0) and X(t f ) unknown. (40) 

Due to the known boundary conditions being specified at the two ends of the 
slewing tiae, this problea is usually called the two-point boundary-value 
problea. This kind of split boundary conditions usually result froa the 
large-angle aaneuver requirements, in which the initial (t=0) and final (t=t f ) 
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states of the system are specified. By solving this problem, we can obtain the 
optimal control (based on the necessary conditions). The often used solution 
strategy is to change the boundary value problem to the initial value problem, 
i.e., find X(0), the missing initial costates. Once X(0) is obtained, one can 
solve the equations (39) as an initial value problem by using any numerical 
integration method. However, owing to the nonlinearity of the equations, there 
is generally no analytical solution to this problem or simple numerical method 
to obtain the solution except for some very simple cases such as the linear 
time-invariant case. The numerical iteration method is the general approach to 
the this problem. 

To start an iteration process, one usually needs an initial guess of 
X (O) (0). Then, equations (39) or their equivalent form (the linearized version 
of (39)) are solved and a x (0> (t f ) is obtained. Based on the difference 
Ax(t f )=x (0) (t f )-x(t f ), the correction to X(0), AX(O), is obtained. This gives 
us a new initial value of X(0), X (1) (0). Hence, the next iteration begins. The 
iteration process can be terminated when |X (k+1) (0)-X <k> (0) | is less than a 
given error limit. One can see immediately that if the beginning guess X (0) (0) 
is close to the true value (converged value) of X(0), the solution will 
converge and less iterations are needed. However, a "good" guess of X(0) is 
often difficult to obtain for the general nonlinear problems. 

Therefore, the effort for solving the TPBVP is two fold. The first is try 
to establish a good iteration (correction) method with a wide convergence 
interval so that it can guarantee convergence even for a "poor" initial guess. 
The other is try to find a "good" initial guess based on the characteristics 
of the practical problem and using some simplified mathematical models. In 
this report, we use the quasilinearization method. He also use the solution of 
X(0) from the simplified linear, time-invariant model of the system as the 
initial guess for starting the iteration process. 

4.1 Linear and Time-Invariant TPBVP 

For linear, time-invariant versions of equations (39) (refs. 1-2), 

z=Az, where z T =[x T , X T ] (41) 
its transition matrix (constant exponential matrix). 
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can be used to obtain the initial costates (closed fora solution) 


X(0)=A;J[x(t f )-A ii x(0)] 


( 42 ) 


4.2 Nonlinear TPBVP 

The continuation (relaxation) process (to increase the participation of 
the nonlinearity in the solution) and the differential correction (for 
determination of the initial costate variables) have been used in references 
1-2 for the 2-D slewing problem. However, as stated in ref. 3, the extension 
of these techniques to the 3-D slewing problem has encountered a numerical 
problem: rank deficiency. 

In refs. 5, 10, and 11, the quasilinearization method has been 
successfully used. In this method, one needs to linearize the differential 
equations ( 39 ) , 


i=g(z), where z T =[x\ X T ], g T =[fJ, f*l (43) 

about an approximate solution of this equation in the following fora (a series 
of linearized, time-variant, nonhomogeneous equations): 

i< k * 1) =Og/a z ) z (k * 1) + h(z (k) ) (44) 

where z (k) is the kth solution of the same linearized equation. It is also the 
kth approximate solution of the original nonlinear equations (43). Here, the 
boundary conditions, (40), are naturally adopted as the boundary conditions of 
the linearized equations, (44). The control expressions, (37), also need to be 
linearized (ref. 8): 

u <k * 1) =u <k) -R- , [P(Aa)] T X (k) -R- , [E+P(« <k) )] T AX (45) 

where Aa=a <k * 1) -a (k) , and AX=X (k4l) -X (k) . By assuming that 
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(46) 


u (k) =-R- 1 [E+F(a (k) )] T X U,) 

for the unbounded control case, equation (45) can be rewritten as, 

u (k*l) = _ R -l [ p( Aa)] T x n«)_R^ [ E + p(a <k, )] T X (k + ,) (47) 

However, in the bounded control case, equations (38) are considered, that is. 


(k). 


-u or u ^ 
ib ib 


u 


(48) 


-{R _ 1 [E+F(a <k) )]V k) } 1 
Accordingly, at the (k+l)st step, u (k+1) can be determined by 


-u or u . if |{R- 1 [E+F(a (k) )] T X (k) } U u 
ib ib 1 1 1 

-{R‘ t [F(Aa)] T X (k) -R" 1 [E+F(a <k, )] T X (k+,) > i 


( k ♦ 1 ) _ J 


u 


(49) 


So far, an iteration process is formed. In each iteration, only a linear TPBVP 
is solved. It is this property that gives this approach the name 
quasilinearization Method. 

The linear TPBVP can be solved by many ready-made methods. One of the 
frequently used algorithms is the method of particular solutions (ref. 9). 
Let m represent the number of the states (also the costates). Equations (44) 
can also be rewritten in the following font, 


i(t)=G(t)x(t)+H(t)X(t) , X(t)=l(t)x(t)+J(t)X(t) 


(50a,b) 


From the theorem of the linear system, any solution equation (50a) can be 
expressed as the linear combination of its mrtl particular solutions, i.e., 

■♦1 . ■** 

x(t)= E 4 c i (t), as long as E c -1 (51) 

i ■ 1 i l»ll 

where c^ are constants and x‘(t) are the ith particular solution vectors. The 
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method begins by integrating equations (50a, b) forward m+1 times, with the 
initial conditions , 


z 1 (0)= 

*(0)1 

1 

0 

0 

, z 2 (0)= 

I 

X 

■* O ^ 

o 

1 

z"(0)= 

x(0) 

0 

0 

, and z >l+t (0)- 

*(o) 

0 

0 

0 


- 0 - 


- o - 


O 

i 


- 0 ■ 


This gives us m+1 particular solutions, x*(t), x 2 (t), • ••> x (t). 

Substituting these solutions into equations (51), and setting t=t , we have 


M4l *+l 

,5, V A V 1 


(52) 


This is a set of m+1 algebra equations for m+1 unknown constants, c^ . By 

assuming the existence of inverse of the coefficient matrix, we can obtain the 

solution, c=[c c • • • c ] T and c . By doing the following manipulation, 

12 a 


z(0)= E c z‘(0)= 
1*1 1 


■ + 1 
E c : 

i-i i 

c 



' *(0) ' 


c 

1 


C 2 

• 

• 

m 

c 

L a 


one immediately realizes that c=\(0), the missing initial costates. 


4.3 Transformation of Attitude States and Costates 

The following procedure is designed to obtain the solution of the 
nonlinear TPBVP. First, the linear TPBVP based on equation (28) is solved and 
a nominal trajectory is produced, in which the control is unbounded and the 
initial costates are calculated by using the transition matrix aethod . Then, a 
converged solution for the linear TPBVP with bounded controls is obtained by 
iterations starting from the previously obtained trajectory. Note that the 
Buler angles are used in all the above computations. 
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Next, to obtain the starting solution for the nonlinear TPBVP , the 3 
Euler angles and the 3 associated costates are transformed to the 4 
quaternions and their 4 costates (from t=0 to t=t f ). Refs. 10-11 provides us 
the following relationship between the quaternions , q(t) and their costates, 
p(t) : 


3 J 


-d -d “d 
12 3 

d “d d 

0 3 2 

d d -d 

3 0 1 

— d d d 
2 i o J 


3 J 


(53) 


where d^ are constants. For the case q(0)=[l 0 0 0] T , we can choose P Q (0) 0* 

d =0. Then, 

o 


[p/0) p 2 (o) P 3 (o)]=[d i d 2 d 3 ]=d T 


(54) 


The vector d can be determined by 


d= 2 [initial Euler angle costate vector] 


(55) 


This result can be proved if one compares the related state and costate 

equations for both linear and nonlinear TPEVPs (for the case w(0)-0) . After 

finding the q(t) by using a nonsingular transformation between q(t) and the 

Euler angles, & (t), 6 (t), and 0 (t), one can use equations (53-55) to obtain 
12 3 

p(t). 

Finally, the nonlinear TPBVP is solved through the quasilinearization 
process and the method of particular solutions [9]. 


5. NUMERICAL RESULTS 

The following parameters of the orbiting SCOLE are used in this paper 

[6]. The inertia matrices of the Shuttle and the reflector, about the mass 

2 

center, o and o , respectively, are (slug-ft ): 
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'905443 

0 

145393 


'4881.375 0 0 

0 

6789100 

0 

» 

0 4969.594 0 

145393 

0 

7086601 


0 0 9921.969 


The material properties of the mast are I EI=GI^=4E+7 lb— ft , pA-0. 09554 

slug/ft , pi =.9089 slug-ft, and L=130 ft. The masses of the Shuttle, the mast, 

and the reflector are (slug) I 6366.46, 12.42, and 12.42, respectively. The 

location of the mass center of the reflector is x —18.75 ft, and y -32.5 ft. 

r r 

The control saturation levels are: | f lx | = | f , y | = l* lz I* 1E+ * ft~lb. 

If | =| f I =| f | =| f Is 10 lb, | f | =1 f I S 800 lb. The first five natural 
frequencies are (hz): .2740493, .3229025, .7487723, 1.244013, 2.051804. 

The numerical tests based on the previously described method have been 
performed for the roll-axis slews, pitch-axis slews, as well as arbitrary-axis 
slews. All these tests are rest-to-rest slews and the iteration process is 
terminated after the initial costates are reached within five digit accuracy. 

The following procedure is designed to obtain the solution of the 
nonlinear TPBVP. First, the linear TPBVP based on Equation (28) is solved and 
a nominal trajectory is produced, in which the control is unbounded and the 
initial costates are calculated by using the transition matrix method. Then, a 
converged solution for the linear TPBVP with bounded controls is obtained by 
iterations starting from the previously obtained trajectory. Note that the 
Euler angles are used in all above computations. Next, the Euler angles and 
the associated costates are transformed to the quaternions and their costates 
[10, 11], to obtain the starting solution for the nonlinear TPBVP. Finally, 
the nonlinear TPBVP is solved through the quasi linearization process and the 
method of particular solutions . 

Case 1 is a 90 deg slew about the roll(x) axis with only three torques on the 

Shuttle as the control, i.e., u=[f f f ] T . The weightings for the state, 

lx ly lz 

0.0. and 0 are chosen to be zero matrices, with the consideration that a 
12 3 

non-zero choice will improve the responses [5]. The control weighting is 
selected as R=Diag[lB-6, IE-6, 1E-6] T , since the small values here imply the 
small costates (which is advantageous for the numerical convergence) for the 
same control (see Equation (37)). The slewing time, t^ =28 sec, makes the slew 
near the minimum-time— slew as used in [5] for the planar SC0LE configuration. 
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The results for this slew are shown in Figs. 3, in which the solid lines 
represent the responses of the linearized system (27), while the dotted lines 
stand for those of the nonlinear system. The three attitude angles are plotted 
in Figs. 3a-c. The roll angle, 0 ^ , for both systems is almost the same, but 
the pitch (y) and yaw (z) angles are different although the magnitudes are 
quite small. Associated with these differences are the differences in the 
controls f and f , shown in Figs. 3j— k. The roll - axis torque, f * (Fig. 

ly lz lx 

3i) is near the bang-bang type. There are little differences in the first 

three modal amplitudes (Figs. 3d-f) between the two systems, but the 4th and 

the 5th modal amplitudes (Figs. 3g-h) present larger relative differences. 

Since the second mode describes mainly the deformation of the mast in the y-z 

plane, which is perpendicular to the slew (x) axis, the second modal amplitude 

has the largest peak value among all the five modes. 

Case 2 is a 90 deg slew about the x axis, but using all 9 controls. <^=0, 

i=l,2,3, and R=Diag[lE-6, IE-6, IE-6, 8E-2, 2E-2, 9E-2, 4E-2, 8E-4, 3E-4] T . 

Due to the increase in the control effort, the slewing time can be shortened. 

t =12 sec is selected numerically by the judgment that the maximum 
f 

displacement of the mast at the reflector end is less than 10% of its total 
length, to be consistent with the small deformation assumption. To make a 
comparison, the slewing results for the rigid spacecraft sxxlel (Equation (27)) 
are also obtained by using the same Q^s and R. Figs. 4 give the results for 
the present slew, where the dashed lines represent the time histories of the 
rigid nonlinear system. 

The three systems have less differences in 0^ (Fig. 4a). and 0^ (Pigs. 
4b-c) are still very small, but the peak values are several times larger than 
those in Figs. 3b-c. The differences between the flexible and rigid nonlinear 
systems are small, but the differences between the flexible nonlinear and the 
the flexible linear systems are relatively large. The reason is that the 
quadratic terms of the angular velocity of the Shuttle, <■>, (Equation (25)), 
have been used both in the rigid and flexible nonlinear systems , but do not 
appear in the flexible linearized system. Therefore, these quadratic terms 
play an important role the 3-D large-angle rapid slewing problems. The similar 
differences in the three systems are also reflected in the control histories 
(Figs. 4i-p). The small differences between the flexible and rigid nonlinear 
systems are caused by the deformation of the flexible mast. For the rapid 


2.16 



slews witH lsr^e control torques nnd forces^ the def 01 us t ion will further 
increase , and so will the differences between these two systems. 

The five modal responses are quite different from those in Case 1. Some 
of the peak values increase about 30 times, and the response phases change. 
For example, before t = t^/2, ot^ is mostly positive in Case 1, but negative in 
Case 2, and this change reverses after t=t y /2. This change is caused by the 
controls at the reflector [5]. 

In the present slew, is the main attitude angle and mode 2 has the 
largest deformation in the y-z plane among all the five modes (Figs. 4d-h) . 
Within Figs. 4, after comparing the responses of the linearized system with 
those of the nonlinear systems, we see that there are less relative 
differences in 0 and a than in 0 (0 ) and other modes. These results imply, 

12 2 3 

for the slew considered here, the major modes (rigid 0 ^ and flexible a z > have 
the largest overall displacements but the smallest relative differences 
between the linear and nonlinear systems. On the other hand, the secondary 
modes (0,0, a, a, a, a) have smaller overall displacements but larger 
relative differences. As a consequence, the major controls (f <y ) (Figs. 4q) 
and all the remaining secondary controls (Figs. 4i-p), as well as the major 
deformation V(z,t) and secondary deformations U(z,t) and d(z,t) (Figs. 4r~t) 
also yield the same results. These results may lead to the following 
conclusion: the linearized equations can represent very well the nonlinear 
equations for the major slewing motion even for large displacements, but not 
so well for the secondary motions, even for small displacements. The 
explanation for this fact might be that the magnitudes of the nonlinear terms 
have a certain level which is less than that of the linear terms 
representative of the major motions, but is great enough to compete with that 
of the linear terms for the secondary motions. 

It should be mentioned that these facts can not be observed in the planar 
slewing problems studied by many authors (for example. Refs. 1, 2, and 5). In 
those researches , the differences between the linearized system and the 
nonlinear system are very small because all the modes are planar modes and, 
hence, the first several modes are all major modes. 


2.17 


Case 3 Figs. 5 show the results for a simultaneous three-axis slew (0^=60, 
0 =30, 0 =45 deg). The weightings for the states are Q 1 = Q 2 = Q 3 =0 * In th * B case » 
The Shuttle torques (f , f , and f ) and the reflector forces (f . and 

f ) are used. The associated weighting for the control is R=DIAG(lE-4, IE-4, 

Ay 

IE-4, 0.6, 1.4-3). The slewing time, t f , is 40 sec. The solid lines m the 
figures 5a-h responses of the rigidized nonlinear system, equation (27), while 
the dotted lines represent the slew results for the flexible nonlinear system. 


6. CONCLUSIONS 

The application of Pontryagin's Maximum Principle to the large angle 
slewing maneuver problem has been extended to the slewing of a 3 D flexible 
spacecraft (SCOLE). A numerical simulation procedure based on the 
quasilinearization algorithm for solving the resulting nonlinear TPBVP has 
been established and tested successfully for several examples. The general 
nonlinear dynamical equations developed here contain all the quadratic terms 
of the angular velocity of the main body and their coupling with the first 
order modal amplitudes and modal rates. It is suggested that higher order 
terms be included if a further analysis is conducted. The numerical results 
show an inportant fact that the linearized system can represent the nonlinear 
system adequately for predicting the major motions but not as well for the 
secondary motions. The quadratic terms (nonlinear) of the angular velocity of 
the main body (Shuttle) cannot be neglected for large-angle rapid maneuvers. 
The differences between the responses of the rigid and flexible nonlinear 
systems are small because the deformation of the flexible part (mast) is 
small. For further research, it is recommended that the applicability of this 
method to more complicated systems be established. 
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Pig. 3. (continued) Case 1, 90 degree roll axis slew. Shuttle torques only. 
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Pig. 3. (continued) Case 1, 90 degree roll axis slew. Shuttle torques only. 
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Fig. 4. Case 2 , 90 degree roll axis slew, all controls. 
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(continued) Case 2, 90 degree roll axis slew, all controls 
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Pig. 4. (continued) Case 2, 90 degree roll axis slew, all controls 
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Pig. 4. (continued) Case 2, 90 degree roll axis slew, all controls. 
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Fig. 5. Results for the 3-D slew. Case 3: 60-30-45 deg-Slew, 
Shuttle Torques + Reflector Forces, 

R=Diag(lE-4, IE-4, IE-4, 0.6, 1.4E-3). 
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Fig. 5. Results for the 3-D slew, Case 3: 60-30-45 deg-Slew, 
Shuttle Torques + Reflector Forces, 

R=Diag(lE-4, IE-4, IE-4, 0.6, 1.4E-3). 
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III. EFFECT OF STRUCTURAL OFFSET AXIAL SHORTENING, AND 

GRAVITATIONAL TORQUE ON THE SLEWING OF A FLEXIBLE SPACECRAFT 


1 . Introduction 

The direct application of Pontryagin's Maximum Principle to the attitude 
maneuvers of spacecraft has been conducted by many authors (Refs. 1-2). 
Recently , some effort has been made to utilize this principle to more 
complicated structures (Refs. 3-4). In Ref. 3, the rapid slewing of a 
2-dimensional flexible orbiting spacecraft, a Shuttle-beam-reflector system, 
has been considered. It is observed (Ref. 3) that the time response history of 
the nonlinear system has a shift from that of the linearized system, but the 
reason for this was not clear. In continuation of this study, the present 
paper will first answer this question by examining the equations of motion and 
by presenting more numerical examples. 

The so called axial shortening effect of a beam induced by its transverse 
displacement has been brought to attention by some authors (Refs. 1-2, 5). 
Although the shortening terms have been included in the equations (Refs. 1-2), 
their effect on the slew lacked quantitative analysis; specifically, the 
numerical examples with and without these terms were not provided. On the 
other hand, a numerical example in Ref. 5 shows that large differences do 
result between models with and without the shortening effect. But the 
numerical example is only for an uncontrolled dynamical response case and the 
main body's motion is prescribed. In the present paper, therefore, the 
shortening terms are considered in the formulation of the equations of motion 
and numerical examples both with and without these terms are presented to 
compare the difference between them. 

Finally, the gravitational torque terms are modeled and included in the 
equations to show their effect on the slewing motion. 
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2. Dynamical Equations 


System Configuration 

The spacecraft model used here (see Fig. 1) is composed of two rigid 
bodies , representing the Shuttle and a reflector (Refs. 3-4, 6), connected by 
a flexible beam through fixed joints, in the plane of the Earth orbit. One end 
of the beam is assumed to be connected to the mass center of the Shuttle, 
while the other end is, in general, not connected to the mass center of the 
reflector. This offset is represented by x in the i direction of the 

A r r 

reflector's coordinate system, (i , k ). It is this offset parameter that will 

r r 

be examined in this paper. 0 is the rotation angle of the Shuttle fixed 

coordinate system, (i , k ), with respect to the orbital coordinate system, 

(i , k ). u(z,t) and <>(z , t )=6u/3z describe the transverse displacement of the 
o o 

beam and the rotation angle of its cross section from its undeformed position, 
respectively. Both u and $ are assumed small and can be expressed by the modal 
superposition formula u=^g i (zjc^ (t) , (z)cc 4 (t) ; where £ is the ith modal 

function, is the ith scaled modal amplitude, and ^^dj^/dz. 

The effect of the offset on the slewing is analyzed by changing the value 
of x^ . Towards this purpose, the partial differential equation for the free 
vibration of this structure has been solved by using the separation of 
variable method. The natural frequencies and modal shape functions have been 
obtained (the assumed mode method was used in Refs. 1-2, 5), for different 
value of x^ . It is observed that the natural frequencies decrease as the 
offset distance increases. 

Kinetic Energy 

The kinetic energy of the system about the mass center of the system, c, 

can be expressed as T=T +(T +T )-T , or 

• b r o 

t = ijj r i 2dm + i[J b+r ,;,2<i ” " 
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where r is the position vector from 0> (Shuttle mass center) to an arbitrary 

mass element in the system and « t is the total mass of the system. The 

integration subscripts. V. "b\ and V. -ean that the corresponding 

integration is throughout the Shuttle, the beam, and the reflector, 

respectively. in Bq. (1) represents the kinetic energy of the Shuttle about 

0 t = -I e 2 ,* where 6 , w is the orbital rate and I # is the moment of 

B 9 b 2 • t t 00 

inertia of the Shuttle, 


Two assumptions for the deformation of the beam have been made: (1) The 
length of the beam does not change; (2) The rotary inertia and the shearing 
force are neglected. Consider now an element dm on the beam (see Fig. 1), 
which has a coordinate, z, before deformation and, z-tz, after deformation 
along the k # axis, where Az is the "shortening" amount due to the deformation 
and can be determined by solving the following functional. 


Z =| Z ^ J l + (§^] dz or approximately, £l+ |[|^) ] dz 


where small (3u/3z) is assumed. Then, 



By dropping the second term, which is a higher order term, we have 


This indicates that Az is a function of z and u(z,t) (a function). Also, 
Az is is a second order term in the modal amplitude, a. Equivalent 
developments for the shortening effect are also presented in Refs. 1 and 5. 
The position and velocity vectors of a mass element dm on the beam and the 
reflector are, respectively, 
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r=ui+(z-Az)k; r=[0 (z-Az)+u]i-(0 u+Az)k; (k, i)-(k , i ); 

t t • • 

r=(u +x)i+(L-AL-x$ )k; and r=[6 (L-AL-x$ )+u ]i-{6 (u +x)+AL+x$ ]k 

r r trrir r 

where L is the length of the beam. T , T , and T in Eq. (1) can be obtained 

b r c 

by using the following equations: 

l b ^ zdn: if o pAz [J o (i?) dy ] dz 

J A Z da=|J p A [f (gj dy]dz 


■HK (*§ M* 2 “ i } 


dz 


After dropping third and higher order terms in modal amplitudes, a, we obtain 

T ) 2 (I+2a T m *hx T M a) )a T (m+M a) *hx t M a 

20 m 2 0 2 4 2 3 

where I, m , m . M , M , and M are constants, vectors, and matrices, 

a 2 2 4 3 

respectively. 


Potential Energy 


The gravitational energy and the elastic energy of the beam are 


V +V = 

• 0 



dz + i“n 

l 0 


k *J'k 

o o 


= 4o T Ko + jju 2 (J sin 2 0 +J cos 2 0 -J sin20) 
2 2 0 11 33 13 


where k =-isifl9 + kcos6, El is the constant flexural rigidity of the cross 

o 

section of the beam, and J is the inertia tensor of the system, with J being 
the functions of the modal amplitudes, a. 
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Generalized Forces 


The shortening effect is also considered in developing the virtual work 
of the controls and the associated generalized forces , and Q^. For example, 
the force arm for a control force is affected by the shortening effect. 

Dynamical Equations 

By using the Lagrangian equations, we can obtain the dynamical equations 
of the system in the following matrix form 


I+2a T m +a T M a j (m +M a) T 


0 

a 2 j 2 4 



m +M a 1 M 


a 

2 4 : 3 


„ m 


-2da T (m^+M 2 a)-a T M 5 a -6V ^/80 +0^ 

-2BM a*6 2 (m +M a)-K a-0V /da +Q 

a a 2 s' a 

where M =M +M . m , M and M are linear functions of z . M„, M, contain the 
4 a 5 a a S r 2 5 

components of the shortening terms. From these equations, we see that the 

terms containing m and M are nonlinear terms of first order in a or a, 

a a 

while the terms containing M , M and M are of second order. Therefore, for 

2 4 5 

moderate nonzero values of x , the influence of the structural offset can be 

r 

greater than the shortening effect. The linearized equations can be obtained 
by neglecting all nonlinear terms. 


( 2 ) 


T • T 

Li"?. 


0 


0 

4* 

~dV JdB 

V 

m ; M 

2 | 3 


a 


-Ka 


-ev Jda 

V 


where "LIN" means constant and linear terms. Note that on the right side of 
Eq. (3), the structural offset and the shortening terms disappear. 
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3. Optical Control 


The optimal control for the slew problem is derived by using the 

techniques similar to those used in Ref. 3 and will not be repeated here. The 

controllers are provided by one control torquer on the Shuttle (u^), one 

control force on the reflector (u^), and two force actuators on the beam (u 2 

at L/3 and u at 2L/3). Each of the controls has its own upper and lower 
3 T 

saturation levels. The cost functional, J=sf u T Ru dt, is used in the present 

** 0 

paper, where u is the control vector, R is the control weighting matrix, and T 
is the slewing time. The resulting two-point boundary-value problem is solved 
by using the quasilinearization algorithm. 

4. Numerical Results 

The parameters of the orbiting SCOLE (Ref. 6) (Spacecraft Control 

Laboratory Experiment) are: EI=4E7 lb“ft 2 , pA=0. 09554 slug/ft, L-130 ft, 

m =6366.46 slug, m =12.42 slug, I =12.42 slug-ft 2 , I =4881.375 slug-ft 2 , 
8 r s r 

g> =0.001 (rad/s) (orbital altitude h*981 km). The natural frequencies (hz) for 
x =0 and x =32.5 (ft) are: 0.3365257, 2.062547, 5.316669; and 0.3199540, 

r r 

1.287843, 4.800169, respectively. All numerical simulations are 90 degree 
rest-to-rest slews and can be represented by: 

Case 1 x =0, u=u , R=1E~6, T=27.6 (s) 

Case 2 x =0. u=fu u u u ] T , R=DIAG(lE-6, .15, .21, IE-4), T=8.196 (s) 
y ’ ^ r 1 2_ 3 4 

Case 3 x =32.5 (ft), u=u^ , R=1E~6 , T=27.6 (s) 

Case 4 x =32.5 (ft), u=[u u u u 1 , R=DIAG(1E“6, .15, .21, 1E~4), 

r l £ J v 

T=8 . 196 (s) 

Pigs. 2a-g display the time histories of 6(t), u(L,t), ♦(L,t), u(t) for 

Case 4. Clearly, the response of 0(t) for both linear and nonlinear systems 

are very close. However, there exist some differences between the two systems 

in u(L,t), 0(L,t) and the controls, u. The difference is primarily due to the 

offset x (here, x =32.5 ft). When x =0, this difference can be reduced 
r r ** 

markedly, regardless of whether the shortening effect and gravitation are 
considered. It is also interesting to know that the controls have large 
differences only around the sud-slew-tism. 
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Table 1 Tip Displacement and Tip Angle 



Max-Disp 

(ft) 

Min-Disp 

(ft) 

Max-Ang 

(deg) 

Min-Ang 

(deg) 

e-Disp 

Case 1 

Linearized 

0.37727 

-0.37727 

0.27402 

-0.27404 

— 


AL=0, (i> =0 

0.37727 

-0.37727 

0.27403 

-0.27404 

0.0% 

x -0, u-u 
r 1 

AL*0, w°=0 

0.37727 

-0.37727 

0.27402 

-0.27404 

0.0% 

T=27 .6 s 

AUO, (*)°/0 
0 

0.37728 

-0.37728 

0.27403 

! -0.27405 

0.0% 

1 

Case 2 

Linearized 

13.072 

-13.072 

9.6216 

-9.6216 

— 


AL=O f u =0 

13.186 

-13.186 

9.7050 

-9.7050 

0.87% 

x ”0 9 u=u 
r 

AL*0, u°=0 

13.154 

-13.154 

9.6847 

-9.6847 

0.63% 

T=8 . 196 s 

AL*0, <0*0 
0 

13.153 

-13.154 

9.6842 

-9.6845 

0.62% 

Case 3 

Linearized 

0.38812 

-0.38812 

0.30342 

-0.30340 

— 


AL=0, <■) =0 

0.38590 

-0.40802 

0.29110 

-0.33940 

5.13% 

x *0, u-u 
r 1 

AL*0, u°=0 

0.38600 

-0.40732 

0.29118 

-0.33884 

4.95% 

T=27 .6 s 

AL*0, u *0 
0 

0.38586 

-0.40803 

0.29094 

-0.33981 

5.13% 

Case 4 

Linearized 

12.191 

-12.191 

9.1082 

-9.1082 

■ ■ 

— 

AL=0, to =0 

12.734 

-12.061 

9.4541 

-9.1299 

4.45% 

x /0, u=u 
r 

AL*0 P u°=0 

12.796 

-12.052 

9.5067 

-9.2030 

4.96% 

T=8 . 196 s 

AL*0, o°*0 
0 

12.795 

-12.054 

9.5052 

-9.2061 

4.95% 


Table 1 lists the maximum (minimum) values for the displacement, u(L,t), 
and angle, $(L,t), of the beam during the associated slews for all cases. The 
first line in each case lists the results for the linearized system, while all 
remaining lines represent those for the nonlinear system with different 
considerations. For example, AL=0 means the shortening effect is not 
considered. The last column gives the largest relative displacement error, 
with respect to the linear results, based on 


(e )* 
Di «p * 


=M.x[|HAX.-M»X u|1 |/|MAX lI11 | . |MIN,-MIH lI1( |/|MIN lI1( |] 


Nonlinear System vs. Linearized System First, let us examine line 1 and line 
2 in each case. In Case 1, since no offset, the differences between the two 
lines are very small. In Case 2, where more controllers are used and the 
slewing time is shortened, the differences increase synmetrically 
(|MAX| = |MIN| ), in spite of x =0. Case 3 uses the same slewing conditions as 
used in Case 1, except X r =32.5 ft. This offset shifts the envelop of the 
response downwards and results in a larger relative displacement error than 
that in Case 2. Case 4 is the combination of Cases 2 and 3. The shift now is 
upwards which is due to the inclusion of more controllers. When more 
controllers are used (Ref. 3), the phase of the response reverses, so do the 
maximum (minimum) amplitudes. 


Shorten Effect By comparing line 2 and line 3 in each case, we can see 
that the shortening terms (1) reduce the amplitude (Cases 2 and 3); (2) 
increase the amplitude (Case 4); and shift the response upwards (Cases 3 and 
4) . These observations coincide with the fact that Az is only a second order 
effect compared with the offset effect. 

Gravitational Effect By observing lines 3 and 4, we can conclude that the 
addition of the gravitational torques into the equations of motion has a very 
small effect on the slew, although they shift the response downwards. This is 
because the orbital rate is much smaller than the slewing rate and the 
magnitude of the gravitational torque term is much smaller than that of the 
active control torque term. 
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5. Conclusion 


Generally , the linearized system can predict the system dynamics very 
well in the slow slewing case. However, in the rapid large-angle slewing 
problem, the responses of the system deviate noticeably from those described 
by the linearized equations if the effects of structural offset and axial 
shortening are included in the simulations. The structural offset (if any) 
results in a first order nonlinear effect. The shortening effect causes only a 
second order nonlinear effect and may not be considered, in the controlled 
simulations, unless the deformation is out of the linear range. The 
gravitational effect can be safely neglected in the slew motions considered 
here . 
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Fig. 1. 2-D Orbiting SCOLE Configuration and Beam Element 
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Pig. 2. Results for the 2-D Slew, Case 4: x=32.5 ft, 
u=[u u u u ] T , t =8.196 sec. 
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Fig. 2. Results tor the 2-D Slew, Case 4: » r =32. 
u=[u u, u t =8.i9S sec. 
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IV. THE EFFECT OF ADDITIONAL DESIGN PARAMETERS ON THE 
LQR BASED DESIGN OF A CONTROL/STRUCTURAL SYSTEM 


I. Introduction 


Some of the difficulties encountered in controlling large space structural sys 
terns arc attributed to their inherent fexibility. As the size of these systems increases, 
due to payload limitations, the total mass cannot be significantly increased. In the 
evolution of such systems, it often becomes necessary to include additional clc 
ments into the design such as those resulting from additional actuators, sensors. 01 

experimental modules 1 . 


Sometimes the resulting control may be optimal only for a prescribed struc- 
tural design. But if we try to later change some of the structural design parameters 
— even by a small amount — then the previously designed control may no longer 
be able to satisfy the mission specifications. If we try to change some of the structur 
al design parameters, perhaps the control system performance will be bctlei than 
before in some sense — i.c. more robust, better transient time constant, reduction 
of initial overshoot amplitudes, etc. On the other hand, a change of some other 
structural system parameters may improve the structural design, but at the expense 
of control system performance. 

In Ref. 2, a combined structural and control optimization problem was for- 
mulated using an optimality criteria approach for the orientation and shape control 
of a free-free beam in orbit. The combined cost function included a form of the 
regulator cost, augmented with the (constrained) weight of the whole structure 
logether with the appropriate Lagrange multiplier. Optimality criteria were derived 
for minimizing the combined cost function and the configuration of the stiucturc 
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obtained was used for the synthesis of control laws using linear quadratic regulator 
theory. The configurations obtained by the combined approach required less con 
irol effort for shape and orientation control of the orbiting beams than that asso- 
ciated with the separately designed structure and LQR-based control systems. 

For the large scale space structural design, it is important to reduce the mass 
of the structural system. The satisfaction of the control requirements during the 
combined control/structural design is also important. The multicritcria optimiza- 
tion approach will be needed for use in the field of structural redesign, which 
should allow a large amount of freedom and variety in selecting the potentially 
large number of design variables. In this sludv, we try to use the quadratic cost 
function and the control properties, such as the transient response time of the 
system attitude, as design criteria. The maximum allowable values of the structural 
mass and the saturation levels of the forces and torques provided by the control 
actuators arc used as the system constraints. A free-free orbiting uniform beam 
with an articulated payload will be considered here as a simple model. The addi- 
tional design of the beam diameter is based on LQR techniques using multicriteria 
which include the cost function and the transient response lime of the attitude 
motion of the beam, subjected to the limited mass of the structure and the satura- 
tion levels of the actuators. The numerical optimization procedure and simulation 
is done using the IBM mainframe computer system. 


II. Equations of Motion 


The dynamics of a long flexible beam in the plane of the assumed circular 
orbit can be expressed (after neglecting the second order effect) as 2 


2 

0 + 3 to 0 

O 



(1) 
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_ t; U 

M Y + K V = FV 

Equation ( 1 ) describes the pilch motion or beam type satellites (rigid body motion) 
and equation (2) gives the vibratory motion of the same beam without damping, 

where 


0 is the pitch angle 
T p is the pitch torque 

0) n is the orbital frequency 

J is the transverse moment of inertia 

\1 is the n x n positive definite mass matrix 

K is the n x n positive definite stiffness matrix 

F is the n x p input distribution matrix 

V is the n x 1 displacement vector 

V is the p x 1 vector of force inputs 


Here the uniform beam with two additional masses in a ciicular oibit is 
considered ( Fig. 1 ). D and d stand for the outer diameter and the inner diameter 

of the beam, respectively. Based on the finite element method \ the beam is as- 
sumed to be divided into four elements with each element having the same geomet- 
rical size and material properties. Five force actuators and one torque actuator are 
assumed to be added to each of the five joints and at joint 3, respectively. The force 
actuators arc assumed to produce forces parallel to the positive Y direction, where- 
as the torque is assumed to act about the Z axis. The two additional masses, m, and 


m 


n ' 


can 


are assumed to be attached to the joints 2 and 4, respectively. The mass matrix 
he expressed by the sum of the beam mass matrix and the attached mass 


matrix, that is 


M = M s + M 


( 3 ) 


The beam mass matrix is represented by ? : 


4.3 


4.4 







m 

0 

f. 



111 R 

0 

M , 





and 

m r = 




- 

0 

0 



0 

0 

Also the stiffness 

matrix of the beam can 

be expressed 


n ) 


(1 ) 
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(k) 
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2/ 
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(k) 
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/ 

-12 -61 

12 

-61 

M 
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_ 61 2 / 2 

-61 4/ 2 


( 6 ) 


The joint between two elements i, j is referred by (ij) and the clement number varies 
as k=1 ,2,3.4; / is the length of each clement of the beam; p is the mass density; and 
E is Young’s modulus. A n and I arc section area and section moment of inertia of 
the beam, respectively, which are dependent on the beam diameters, D and d. 


Equations (I) and (2) can be written in slate variable form as: 

X = AX + B U (?) 
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where the 2(10+l)x1 vectors. X. X. anti the 6x1 vector. It. are given hy: 

x = [ e y 1 e v ' ]' 

x = [e v' o v 

I' = [T,, V, V, v, V, vj' 

The stale matrix. A. and the control matrix, B. arc given by 


0 

I 

(Ixll) 

(11x11) 

-3( oj 0 


- 1 

0 

0 -M K 

(1 Ixl 11 


0 

(11x6) 

1 

J 0 

0 M F _ 


(9) 


HI. Additional Design Based on LQR 


Based on LQR theory \ in general, the quadratic cost function is defined by 


r, •• • y •'>/ IS 

OF rDOtf QuAUTY 
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.1 


( 10 ) 


r w i 

= I ( X Q X + 

' n 


U R U) ell 


subject In the state equation of the system from equation (7) 

X = AX + B V 

The cost function, J, is minimized under the optimal feedback control gi\en by 


where P is the positive definite solution of the steady-stale Riccali equation, which 
is 

PA + a'p- PBR b'p + Q = 0 ( 1 21 

From equations (4) and (5). we find the control system properties not only 
depend on the weighting matrices, Q and R, but also depend on the parameters of 
the structure and the actuator locations. If the structural parameters are fixed, the 
regulative range may be very limited for control design. It is possible that the 
change of structural design parameters could extend the regulative range and make 
the control and structural design satisfv the mission requirements. 

Now we consider a class of additional design parameters, Z. which could be 
varied while we design the control system. This means the control, U, is now a 
function of Q, R, Z and t. So the cost function. J, also is a function of Q, R and Z 

Using the extremum principle, assuming that Z results in an extremum of the 
system, we have as a necessary condition for the minimization of J 


dJ ( Q R. Z) 

bZ 


Z 


= 0 


* 


= z 


( 13 ) 
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subject to some specin! control properties 

g ( X. Z. t ) = g n (H) 

where g n is given from the design (mission) requirements. 

From equation (13). Z is obtained, for special eases analytically, or, more 

+ 

generally, by numerical means. The extremum point. Z , may be not unique; thus, 
all extremum points must be compared in order to find the point at which the 
structural performance is optimal, such as the minimization of the total mass of the 
structure. Then, the optimal value of the structural design parameters may be in- 
corporated into the control design to obtain the optimal control. 

In the given simple model (Fig. 1), for example, we may select the section 
diameters, D and d. as the redesigned system variables and fix the transient re- 
sponse time of the pilch motion or the beam. Using the feasible directions of the 
search approach, the optimal values for the diameters would be obtained for the 
given system model and other parameters. For a different given model or parame- 
ter values, the optimization solution would be expected to yield different results. 


IV. Design Multicriteria 


In the LQR process, the assignment of Q and R values normally results from 
a step by step numerical search procedure. For actual LSS design, the interesting 
properties arc the required mass of the structure, the system response time for 
damping rigid motions, and suppressing vibration, as well as the maximum values 
for the actuator forces, etc. A knowledge of the possible expected range of the 
boundary conditions will be helpful in selecting the weighting matrices. Q and R. 
For the combined control/structural optimal design, some design variables for both 
the structure and actuators should be considered. According to these, the control/ 
structural design multicriteria for the system with additional design parameters, Z, 
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could he expressed by 


r r 

I ( Q. R, Z ) = | X ( Z. () Q X ( Z, I ) + 


o 

i 


U ( Z, t) R U ( Z, t ) I dt 

(1?) 

d.J ( Q, R, Z) 

= o 

dz 

z = z 

(16) 

subject to 


X( Z t ) = \ ( Z ) X( Z, t ) + B ( Z ) U( Z, l ) 

(17) 

g ( X. Z, 1 ) = g n 

(18) 

M( Z ) < M max 

(19) 

|U Z. t) | < Umax 

(20) 


where g n is a design requirement 

Mmnx is the limited mass of the structure 

Um, ix is the vector of maximum actuator force 

Z is a vector of time-invarient design parameters 

When g n , Vf max, Umax and the structural configuration arc given, we can 

+ * 

design the optimal structure, Z . and optimal control, U . using the additional 
design formulation based on LQR techniques 


; ' 5 $ 
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V. Numerical Analysis 


We have analyzed and designed an orbiting free-free beam using this tech- 
nique. with pitch and other in-plane flexible modes included as degrees of free- 
dom (see Fig. I). Here vve considered two cases. Case 1 is a uniform solid beam 
(inner diameter d=0). Case 2 is a uniform tubular beam whose wall thickness is 

assumed constant (l)-d=0.04 ft). The material density of 200 Ib/ft 3 . Young's 
modulus of 6E+0 lb/f( 3 . length of the beam of 130 ft and six actuator locations 
arc assumed equal lor both cases. Also the initial pitch angle. 0 o . is assumed 6 

degrees (0.105 rad) and the other initial state components and rates arc assumed 
zero. Without losing general meaning, the weighting matrix. Q. is assumed con- 
stant and a diagonal matrix, that is 


Q = trace I 1000 1000. 100 100 ] 

i i i 1 1 


and the weighting matrix. R. is assumed a unitary matrix multiplied by a variable 
coefficient. By regulating the weighting matrix. R, the response time for the rota- 
lional motion of the beam will satisfy the design requirement (here assumed ten 
seconds). Then the cost function values which depend on the diameter of the beam 
can be calculated by the multicriteria given by equations ( 1 5)— (20) 


First, we suppose the additional payload mass to be 100 1b and to be added 
symmetrically to the beam at joints 2 and 4. The variation of the cost function with 
the outer diameter for the two cases is shown in Figs. 2 and 3. The first cost 
function extremum points for case 1 and case 2 arc deteimined as D-0.45 ft and 
D=0. 32 ft. respectively. If the beam diameters at these points satisfy the mission 
requirements, the structural design is optimal and so is the corresponding LQR 
based control design. Otherwise, the structural configuration, actuator locations or 
material properties should be changed. Comparing the two cases, we find the extre- 
mum points to be different, such that the op-timal diameter for the solid beam is 
larger than that of the tubular beam. In this situation, the solid beam is heavier than 
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the tubular beam for the same required response lime. According to the minimum 
mass requirement for large space structural design, the tubular beam may be much 
better because its mass is much less than that of the solid beam. From Fig. 2, ive 
also can imply that the system may contain an additional extremum point for the 
cost function and its value may be less than the first extremum. If this happens, the 
first extremum point value of the diameter of (he beam should be selected as the 
optimal solution. 

Second, we try to increase additional payload mass up to 1,000 lb and still 
maintain symmetry. There also exist extremum points (see Figs. 4 and 5). but the 
diameter values corresponding to the extremum points differ from those when the 
pavload mass is 100 lb. Comparing Figs. 2 and 4, or Figs. 3 and 4, we find the 
optimal diameter of the beam increases when the additional payload mass in- 
creases. This result indicates that the optimal values of the structural parameters arc 
dependent on the payload added. Table I lists the optimal diameter values and 
their payload ratio for the uniform solid and tubular beams with three different 
payloads: 100 lb. 200 lb and 1,000 lb. It is clear that the payload ratio of the 
tubular beam is greater than that of the solid beam for the same control require- 
ment (ten seconds of the transient response time). The payload ratio of the solid 
beam decreases when the payload increases, but the payload ratio of the tubular 
beam increases. 

When additional payload masses are added asymmetrically with respect to 
the center of the beam, there are no big differences between the cost functions for 
(he symmetric and asymmetric payloads (see Table 2). This may be explained by 
the fact that the incremental moment of inertia about the Z axis due to the payload 
masses is designed to be the same for both the symmetrical and asymmetrical distri- 
bution of the additional payload. Thus for practical system design it is probably 
useful to emphasize the symmetrical distribution of the additional payload wherev- 
er possible. 

Figs. 6-9 provide the transient responses of pitch angle, two deformations, 
and corresponding control torque and forces for the nearly optimized uniform 
tubular beam (Fig. 3). Here the outer diameter of the beam is 0.3 ft which is near 
the optimal value The wall thickness is still maintained at 0.02 ft. The transient 
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response time of pitch angle and all deformations are about ten seconds. Since the 
additional payload masses on the beam arc symmetrically distributed, the deforma- 
tions at joints 4 and 5 arc the same as the deformations at joints 2 and I , respective- 
ly, but their directions are opposite to each other. The maximum torquerequired is 
less than 5,000 ft— lb. All actuator forces are small and do not exceed 25 lb. If the 
solid beam is used instead of the tubular beam, for optimization, the control system 
requires 6,905 ft— lb maximum torque and 133 lb maximum force, which arc 
larger than those for the tubular beam. 


VI. Conclusion 


This paper reviews the effect of additional design parameters on the LQR 
based optimal design of space structural system. A multi-objective cost function 
which includes a form of the standard LQR regulator cost and its partial variation 
with respect to the additional design parameters is considered. The constraints are 
extended to the desired control properties. The optimal multicriteria are derived by 
minimizing the cost function and setting the variation of the cost function with 
respect to (he design variables to zero. This approach is used to determine the 
optimum diameter of an orbiting free-free uniform beam with additional payload 
masses added. From the numerical results for two design models — a uniform solid 
beam and a uniform tubular beam, with two typical additional payloads added 
symmetrically and asymmetrically about the center of the beam, it is found that the 
optimal diameter occurs at the first extremum point of the variation of the cost 
function with respect to the diameter. It is also found that (he tubular beam is 
superior to the solid beam for meeting both minimum mass requirements as well as 
desired transient and control requirements. The study proves that the multicriteria 
design approach should give better results from both the structural designer’s and 
the control designer’s points of view. 
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Table I: Comparison of Optimal Diameter and Payload Ratio of Uniform Solid 
and Tubular Beams for Payloads: 100, 200 and 1.000 lb 


Payload 

(lb) 

Diameter (ft) 
Solid Tubular 

Payload Ratio 
Solid Tubular 

100 

0.45 

0.32 

2.4 C 7 

22 r T 

200 

0.75 

0.42 

2.0^7 

30<T 

1 ,000 

1.60 

1.50 

1.9<* 

30<~7 

Table 2: Comparison of Cost Function Value Varied with Diameter, D, for Sym- 
metric and Asymmetric Payloads (500 & 500 lb; 900 & 100 lb) for the Uniform 
Solid and Tubular Beams 

Diameter 

Solid 

Beam 

Tubular Beam 

D 

Jsvm 

Jasym 

Jsym 

Jasym 

1.3 

30.025 

30.024 

34.437 

34.530 

1.4 

28.902 

28.901 

33.939 

33.937 

1 .5 

27.525 

27.524 

33.845 

33.843 

1.6 

26.505 

26.504 

33.856 

33.853 

1 .7 

27.942 

27.941 

33.865 

33.863 

1.8 

29. 152 

29.151 

33.875 

33.873 

1 .9 

30.778 

30.777 

33.885 

33.883 
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Fig. 1 


Configuration of Orbiting Free-free Uniform Beam with two Additional Payload Masses 
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Fig. 2 Cost Function vs Diameter of the Uniform Solid Beam (here d=0) when two 
50 lb Payloads are Added to the Beam at Joints 2 and 4, Symmetrically 
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Fig. 3 Cost Function vs Diameter of the Uniform Tubular Beam (here D-d=0.04 
ft) when two 50 lb Payloads are Added to the Beam at Joints 2 and 4, Symmetrical ly 
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Fig. 5 Cost Function vs Diameter of the Uniform Tubular Beam (here D-d=0.04 
ft) when two 500 lb Payloads are Added to the Beam at Joints 2 and 4, Symmetrieally 
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Fig. 7 Response of Deformations Parallel to the Y Axis at Joints 1 and 2 for (he 
Uniform Tubular Beam with two 50 lb Payloads at Joints 2 and 4, Symmetrically 
(here D=0.3 ft and D-d=0.04 ft) 
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Fig. 8 Control Torque Trajectory at Joint 3 in Z direction for (he Uniform Tubular 
Beam with two 50 lb Payloads at Joints 2 and 4, Symmetrically [here D=0.3 ft and 
D-d=0.04 ft] 
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Fig. 9 Control Force Trajectories at Joints 1 and 2 Parallel to the Y Axis for l he 
Uniform Tubular Beam with two 50 lb Payloads at Joints 2 and 4, Symmetrically 
[here D=0.3 ft and D-d=0.04 ftj 
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CONCLUSIONS AND RECOMMENDATIONS 


Pcntryagin ' s maximum principle has been applied to study the 
rapid slewing of the SCOLE orbital configuration m three 
dimensions . The quasi linearisation technique for solving the 
resulting nonlinear two point boundary value problem has been 
successfully used for several different examples. The results 
indicate that the linearized system model can represent the 
nonlinear system adequately for simulating the major motions , but 
not as well for the secondary motions. The nonlinear quadratic 
terms of the main body (Shuttle) angular velocity can not be 
neglected for large-amplitude rapid maneuvers. The differences 
between the rigid and flexible nonlinear system responses are 
small because the flexible vibrations are successfully suppressed 
during the maneuvers simulated here. 

In addition it is seen that the structural offset of the 
SCOLE mast attachment to the reflector is associated with a first 
order nonlinear effect. The mast shortening is associated with 
only a second srder nonlinear effect and should be considered 
when the mast deformations are outside of the lienar range. 
Gravitational-gradient effects may be safely neglected for all 
rapid slew maneuvers considered. 

For further research it is recommended to extend the 
application of the maximum principle and the two point boundary 
value problem to mere complicated systems proposed for the future 
CSX program. 

A multi-objective cost function which includes a form of the 
standard LQR regulator cost and its partial variation with 
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respect to additional design parameters is studied here. This 
approach can be employed to determine the optimal diameter of a 
free-free orbiting uniform beam with additional payload masses 
added when constraints are placed on the maximum total mass, 
control saturation levels, and transient settling times. This 
study proves that the muiticnteria design approach should 
produce superior results as compared with combinations of 
separate structural and control system design approaches. 
Extensions are recommended to consider more complex systems 
representative of proposed candidate CSI systems. 

The current < 1 9 9 0 — ? 1 ) grant work has been redirected to lend 
greater support to the Controls /Structures Interaction (CSI) 
program and focuses on specific CSI evolutionary configurations. 
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